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ABSTRACT 

We report the discovery of a probable CH star in the core of the Galactic globular 
cluster M14 (= NGC 6402 = G1735-032), identified from an integrated-light spectrum 
of the cluster obtained with the MOS spectrograph on the Canada-France-Hawaii 
telescope. From a high-resolution echelle spectrum of the same star obtained with the 
Hydra fiber positioner and bench spectrograph on the WIYN telescope, we measure 
a radial velocity of —53.0 ± 1.2 km s“^. Although this velocity is inconsistent with 
published estimates of the systemic radial velocity of M14 (e.g., Vr ~ —123 km s“^), we 
use high-precision Hydra velocities for 20 stars in the central 2'.6 of M14 to calculate 
improved values for the cluster mean velocity and one-dimensional velocity dispersion; 
—59.5 ± 1.9 km s“^ and 8.2 ± 1.4 km s“^, respectively. Both the star’s location near 
the tip of the red giant branch in the cluster color magnitude diagram and its radial 
velocity therefore argue for membership in M14. Since the intermediate-resolution 
MOS spectrum shows not only enhanced CH absorption but also strong Swan bands 
of C 2 , M14 joins to Cen as the only globular clusters known to contain “classical” 

CH stars. Although evidence for its duplicity must await additional radial velocity 
measurements, the CH star in M14 is probably, like all field CH stars, a spectroscopic 
binary with a degenerate (white dwarf) secondary. The candidate and confirmed CH 
stars in M14 and lv Cen, and in a number of Galactic dSph galaxies, may then owe 
their existence to the long timescales for the shrinking and coalescence of hard binaries 
in low-concentration environments. 


Subject headings: stars: carbon — stars: binaries: — Galaxy; stellar content — 
globular clusters; individual (M14) — galaxies: Local Group — galaxies: stellar 
content 


1. Introduction 

Among the presently known types of carbon stars, only the CH stars have abundances and 
kinematics which are indicative of membership in the Galactic halo (McClure 1985). Since stars 
of mass M ~ O. 8 M 0 are not thought to experience the third dredge-up mechanism during their 
ascent of the asymptotic giant branch (e.g., Iben 1975), the origin of the enhanced abundances 
of carbon and s-process elements in these giants remained a puzzle until the discovery that all 
field CH stars are binaries composed of a red-giant primary and a degenerate (i.e., white dwarf) 
secondary (McClure 1984; McClure & Woodsworth 1990). The peculiar abundances of these 
objects are therefore thought to be the result of mass exchange via stellar winds or Roche-Lobe 
overflow during the ascent of the white-dwarf progenitor up the asymptotic giant branch (Han et 
al. 1995). Clearly, the discovery of globular cluster CH stars would have important implications 
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not only for nucleogenesis and globular cluster abundance anomalies (Pilachowski et al. 1996), 
but also for the formation, evolution and destruction of binaries in dense environments (Hut et al. 
1992). 

However, searches for CH stars in several globular clusters based on spectroscopic surveys 
of brightest cluster members (Harding 1962), direct imaging through intermediate-band filters 
(Palmer 1980; Palmer & Wing 1982), and transmission-grating slitless spectroscopy (Bond 1975), 
have proved, by and large, unsuccessful. At present, a handful of stars having enhanced carbon 
and s-process elements have been reported in each of u) Cen (e.g., Harding 1962; Dickens 1972; 
Bond 1975), M22 (Hesser, Hartwick & McClure 1977; McClure & Norris 1977; Hesser & Harris 
1979), M55 (Smith & Norris 1982) and M2 (Zinn 1981). However, while the spectra of these 
stars are characterized by abnormally high CH absorption compared to other cluster giants, they 
usually do not show strong Swan bands of C 2 , suggesting that their anomalous carbon abundances 
probably arise through a different mechanism, such as incomplete CN processing (Vanture & 
Wallerstein 1992), than that operating in “classical” CH stars. Indeed, among this sample of 
CH-enhanced stars in globular clusters, only two are likely to be genuine CH stars. Both of these 
stars, RGO 55 (Harding 1962) and RGO 70 (Dickens 1972), are found in u) Cen. 

In this Letter, we report the serendipitous discovery of a carbon star in the core of the 
poorly-studied globular cluster M14 and argue that it is likely to be a “classical” CH star: a post 
mass-transfer binary consisting of a red giant primary and a white dwarf secondary (McClure & 
Woodsworth 1990). 


2. Observations 
2.1. MOS Spectroscopy 

During an observing run in May 1996 intended to measure radial velocities and Mg b line 
strengths for globular clusters surrounding the supergiant elliptical galaxy M87, we used the 
MOS imaging spectrograph on the Canada-France-Hawaii telescope to obtain long-slit comparison 
spectra of the integrated light of several Galactic globular clusters. During the integrations, the 
telescope was allowed to drift, typically through a spatial extent of about one arcminute, in a 
direction (S ^ N) perpendicular to the slit in order to generate a spectrum representative of the 
integrated light of the cluster. All spectra were accumulated using the B400 grism and STIS2 
CCD which combine to give a dispersion of 3.6 A pixel“^ and a resolution of ~ 8 A. A blocking 
filter was used to restrict the wavelength region to 4500 - 5700 A. 

Unexpectedly, the long-slit spectrum of the integrated light of the Galactic globular cluster 
M14, obtained on 19 May 1996, showed strong Swan bands of C 2 indicating the presence of a 
carbon star in the cluster core. Because of the driftscaning technique employed, the immediate 
identification of the carbon star was in some doubt. By good fortune, however, the carbon star 
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has few bright neighbors to the north or south, allowing a probable identification of the candidate 
based on its E-W position (see Figure 1). On 20 May 1996, we obtained a 600s spectroscopic 
exposure of the suspected carbon star, this time using a l"xl0'^ slit. The blocking filter was 
removed to extend the spectral coverage to the region 4000 - 8000 A. After standard preprocessing, 
the carbon star spectrum was traced, extracted and wavelength calibrated using IRAF.® The final 
spectrum is shown in Figure 2. 


2.2. MOS Photometry 

Of the three published color-magnitude diagrams (CMDs) for M14 (e.g., Smith-Kogon, 
Wehlau & Demers 1974; Shara et al. 1986; Margon et al. 1991), only the photographic study 
of Smith-Kogon, Wehlau & Demers (1974) includes stars brighter than the cluster horizontal 
branch. Unfortunately, these authors did not publish magnitudes, colors or finding charts for their 
program stars. Therefore, to determine the location of the carbon star in the cluster CMD, we 
obtained a number of B and V images of M14 using MOS on 19 and 22 May 1996 (Figure 1). All 
frames were then reduced using DAOPHOT II (Stetson et al. 1990). Empirically we found that 
the CMD based on the two poorest seeing frames (FWHM ~ l".l) exhibited less scatter, reflecting 
the severe MOS undersampling (C/.^44 pixels) in the other frames. Our instrumental magnitudes 
were then calibrated using eight unsaturated local photoelectric standard stars (Smith-Kogon, 
Wehlau & Demers 1974). These stars span a magnitude range of 12.46 < V < 14.59 and a color 
range of 1.17 < B — V < 2.05. The internal random errors in the calibration are ~ 0.02 in V and 
~ 0.04 in (B — V). 


2.3. Hydra Spectroscopy 

Although the carbon star’s small offset of 13^8 ~ 0.3rc from the cluster center (Shawl & 
White 1986; Trager, King & Djorgovski 1995) and its position near the tip of the cluster red giant 
branch (see Figure 3) suggest that it is physically associated with the cluster, radial velocities for 
both the carbon star and M14 are needed to provide unambiguous evidence of membership.® In 
particular, radial velocities are needed to eliminate the possibility that the object is a carbon-rich 
dwarf residing in the solar neighborhood (e.g.. Green et al. 1994), although the a priori probability 
of this is low: to a limiting magnitude of U = 18, the surface density of faint, high-latitude carbon 


®IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under contract to the National Science Foundation. 

®For instance, RGO 153, a carbon star which lies only 18', or (Trager, King & Djorgovski 1995), from the 
center of u Cen, was shown by Smith & Wing (1973) to have a radial velocity which differs from that of the cluster 
by more than 250 km s“^. For comparison, the confirmed CH stars, RGO 55 and RGO 70, are located 9' and 19', 
respectively, from the cluster center. 
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stars (of which dwarf carbon stars make up ~ 13%) is ~ 0.02 deg“^ (Green et al. 1994). Moreover, 
the MOS spectrum (Figure 2) shows only a weak C 2 bandhead at 6191 A, a feature which Green 
et al. (1992) suggest to be unusually strong in the spectra of dwarf carbon stars. Finally, if the 
star in question were a nearby dwarf, it might be expected to show a measureable proper motion. 
Gomparison of our CFHT images with a photographic plate taken in 1952 with the Mt. Wilson 
telescope (see Figure 1 of Wehlau & Froelich 1994) shows no evidence for such a proper motion. 

On 18-19 July 1996, we used the Hydra fiber positioner and bench-mounted echelle 
spectrograph on the WIYN telescope to measure 505 radial velocities (median precision ~ 1.3 
km s“^) for 493 stars in a ~ 0.5 deg^ field centered on M14. In the course of this survey (Cote 
& Welch 1997), we also obtained an 1800s spectroscopic exposure of the M14 carbon star at a 
resolution ~ 0.33 A over the spectral range 4983 - 5252 A. After the program spectra were traced, 
extracted and wavelength-calibrated using IRAF, they were cross-correlated against a template 
spectrum of HD223094, a K5III radial velocity standard. In the case of the carbon star spectrum, 
the adopted template was HD156074, an R star which is known to have a constant heliocentric 
radial velocity of —12.99 km s“^ (McClure 1996). 


3. Discussion 
3.1. Membership 

From our Hydra spectrum of the carbon star, we measure a radial velocity of Vr = —53.0 ± 1.2 
km s“^ (heliocentric Julian date = 2450282.761). How does this compare to the systemic velocity 
of M14? Available radial velocity measurements for this cluster have been reviewed by Webbink 
(1981), who quotes a weighted mean velocity of tv = —123 ± 5 km s“^. This value is based 
primarily on low-dispersion spectra of the integrated cluster light obtained by Mayall (1946), 
as well as a small number of radial velocities (spanning the range -|-9 to —153 km s“^) for four 
Type H Cepheids in the field of M14 accumulated by Joy (1949). An image-tube spectrogram 
obtained by Hesser, Shawl & Meyer (1986) of the integrated cluster light yielded a radial velocity 
of Vr = —25 ± 14 km s“^. Given the obvious difficulties with field star contamination in this 
direction (Z^^ = 21.3°, = 14.8°), as well as the large uncertainties of the Mayall (1946) velocities 

(typical internal error ~ 33 km s“^, but potentially much larger external errors; see Hesser, Shawl 
& Meyer 1986), we have used our new sample of Hydra velocities to determine improved estimates 
of the systemic velocity and internal dispersion of M14.^® 


^°The result of Hesser, Shawl & Meyer (1986) is most likely due to contamination of their spectrum by field stars, 
since the distribution of radial velocities in this crowded field shows a broad peak near ~ —20 km s“^ (Cote & Welch 
1997). Similarly, if the observations of Joy (1949) were made while the Type II Cepheids were near maximum light, 
then the mean velocity of —126 km s“^ he obtained may be (partly) explained by atmospheric pulsations, which can 
approach 30-40 km s“^ in stars of this type. 
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We define a sample of probable M14 members by restricting ourselves to those stars which 
lie within two half-light radii of the cluster center (i.e., 2rh — 2l6; Trager, King & Djorgovski 
1995). This sample consists of 20 stars which have a median velocity uncertainty of 1.2 km s“^ 
and span the range —72.0 < Vr < —47.6 km s“^. BV photometry for these stars (see Figure 
3) demonstrates that they are all likely to be cluster members since they define a smooth red 
giant branch extending from V = 14.38 to 16.23. From this sample, we find a mean velocity of 
Vr = —59.5 ± 1.9 km s“^ and a one dimensional cluster velocity dispersion of Uc = 8.2 ± 1.4 km s“^ 
using the maximum-likelihood estimators of Pryor & Meylan (1993). Based on the close agreement 
between the radial velocity of the carbon star and that of M14, we conclude that the star is indeed 
a cluster member. We also note that, if the star in question is similar to the field CH stars studied 
by McClure & Woodsworth (1990), an additional velocity component due to the orbital motion of 
the red giant primary around the center of mass of the system is expected; among the eight CH 
stars monitored by McClure & Woodsworth (1990), the velocity semi-amplitudes ranged from 4.3 
to 12.1 km s“^, with a mean of 8.1 km s“^. 

3.2. Implications 

Radial velocity surveys of field CH stars have provided compelling evidence that all of these 
systems are post mass-transfer binaries. We therefore conclude that the carbon star in M14 is 
very likely to be a spectroscopic binary having a white dwarf secondary, although confirmation of 
duplicity must await additional radial velocity measurements.^^ It is perhaps notable that, as the 
only two clusters known to contain candidate or confirmed CH stars, M14 and oo Cen are both 
massive, low-concentration systems. A possible connection between cluster concentration and 
CH enhancement was noted previously by McClure & Norris (1977) who, in a prescient remark, 
suggested that “searches for CH stars in the low-concentration clusters M14 and NGC 2419 might 
be profitable”. 

Does environment play a role in the evolution of globular cluster CH stars? The dominant 
dynamical processes affecting binaries in globular clusters are the disruption of “soft” binaries 
through stellar encounters and the shrinking of the orbits of “hard” binaries via energy exchanges 
with intruder stars (e.g.. Hut et al. 1992). Equation (1) of Hills (1984) can be used to estimate the 
separation, Ug, of the widest cluster binaries which are expected to have escaped disruption over 
a Hubble time. Since most Galactic globular clusters have three-dimensional velocity dispersions 
cTy ^ 15km s“^ (Pryor & Meylan 1993), we find ^ 6 AU. This is comparable to the separation 
of the widest binaries in the McClure & Woodsworth (1990) sample of field CH stars (i.e., 1 ^ a ^ 
4.5 AU). Therefore, we conclude that the process of CH-star disruption is unlikely to be important 
for most clusters. 

^^This connection between the field and globular cluster CH stars has been further strengthened by the recent 
demonstration that both RGO 55 and RGO 70 in u) Cen are spectroscopic binaries (Mayor et al. 1996). 
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On the other hand, the shrinking of hard binaries may play a more significant role. At the 
cluster center, intruder stars will shrink the orbits of hard binaries at the rate 


d Ina —2TrGapc 


dt Oc 

where a is the semi-major axis of the binary, is the central one-dimensional velocity dispersion 
and Pc is the central mass density (Hills 1984; Phinney 1996). Integration of this equation yields 
the initial semimajor axis of the binary whose size is halved over the lifetime, to, of a cluster: 


ah = 


27rGpcto' 


( 2 ) 


The orbits of hard binaries that are initially wider than this shrink rapidly, and have a ~ o/i after 
a Hubble time. Conversely, tighter binaries evolve on rather longer timescales, and are largely 
unaffected. Thus, ah is the characteristic hnal separation of surviving cluster binaries (Phinney 
1996). It is worth bearing in mind, however, that equation (2) refers to the cluster center; the 
situation for the cluster as a whole is undoubtedly more complicated. 


In Figure 4 we plot the distribution of at against absolute magnitude, My, for all Galactic 
globular clusters (filled circles) having measured velocity dispersions and central mass densities 
(e.g., see Pryor & Meylan 1993). Absolute magnitudes are taken from Djorgovski (1993). M14 and 
uj Cen are shown as open stars. The open squares show the location of the Galactic dSph galaxies 
(Irwin & Hatzidimitriou 1995). At least some of the carbon stars seen in several of these galaxies 
(i.e., Ursa Minor, Draco, Sculptor and Garina) are likely to be CH stars (McGlure 1984; Aaronson 
& Olszewski 1987; Armandroff, Olszewski &: Pryor 1995). The dashed line shows the semi-major 
axis of the shortest-period binary in the McGlure & Woodsworth (1990) sample of field GH stars 
(i.e., a ~ 1 AU). This size is similar to the mean radii of low-mass, thermally-pulsating AGB 
stars (e.g., see Figure 2 of Boothroyd & Sackman 1988). Closer initial binary orbits would disrupt 
the normal evolution of the carbon donor before it reached the thermally-pulsating AGB phase, 
preventing the dredge-up of carbon. Therefore, Figure 4 suggests that the low concentrations of 
M14 and iv Cen correspond to small values of Pcjac and hence to ah large enough to accomodate 
CH systems with separations typical of field CH stars. 

The discovery of a CH star in M14 therefore provides additional support for the notion 
(McClure 1984) that the process of CH star formation depends on environment, in particular, 
on the rate of hard binary shrinking through stellar encounters. Given the incomplete and 
inhomogenous nature of existing surveys for CH stars in globular clusters, a renewed effort to 
discover such objects in an expanded sample of clusters (particularly those objects which he above 
the dashed line in Figure 4) may provide new insights into the influence of environment on binary 
evolution. 


The authors thank the staffs of the CFH and WIYN telescopes for their outstanding support. 
Thanks also to Robert McClure and Amelia Wehlau for helpful discussions. The research of DAH 
and GLHH is supported through operating grants from NSERC. 
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Fig. 1.— B image of M14 taken with MOS on the CFHT. The carbon star is identified by the small 
circle, while the larger circle denotes the cluster core radius (rc = 50'' according to Trager, King & 
Djorgovski 1995). The cross shows the location of the cluster center (Shawl & White 1986). North 
is up and east is to the left on this image, which measures 3f75x3'5. 

Fig. 2.— MOS spectrum of the newly discovered M14 carbon star. The spectrum has a dispersion 
of 3.6 A pixel“^ and a resolution of ~ 8 A. A number of the most prominent spectral features are 
indicated. 

Fig. 3.— Color-magnitude diagram (CMD) for M14 based on a pair of BV images taken with MOS. 
Both images have FWHM ~ I'.'l, corresponding to 2.5 pixels. The location of the carbon star in 
the CMD (i.e., V = 14.85, B — V = 2.50) is indicated by the open square. Those stars having 
radial velocities measured with Hydra and lying within 2^6 of the cluster center are indicated by 
the filled circles. 

Fig. 4.— Maximum separation of hard binaries, ah, due to shrinking via energy exchanges with 
intruders plotted against absolute magnitude. My, for Galactic globular clusters (filled circles) and 
dSph galaxies (open squares). The dashed line shows the semi-major axis of the shortest-period 
binary in the McClure & Woodsworth (1990) sample of field CH stars (i.e., a ~ 1 AU). M14 and 
LO Cen are shown as open stars. Based on its location in this diagram, NGC 2419 is a prime 
candidate for a search for CH stars. 
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